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Abstract 

There are several solutions to the Einstein equation exhibiting naked-singularity for- 
mation in gravitational collapse, which could possibly serve as counterexamples to the 
cosmic censor hypothesis. It has not been examined seriously, however, whether or not 
such naked singularities are indeed visible to asymptotic observers in astrophysically re- 
alistic stellar collapse, such as the core collapse and the delayed collapse in the final stage 
of its evolution. In this paper, we set the spherically symmetric mass distribution which 
is motivated by the astrophysical scenarios and evolve it with vanishing pressure for sim- 
plicity. We show that even if a naked singularity could appear at the centre, the naked 
singularity is hidden behind an event horizon and cannot be observed by a distant ob- 
server. With this illustration, we argue the existence of an "astrophysical censor" that 
prohibits singularities formed in astrophysical gravitational collapse from being observed 
by a distant observer. 



Contents 

1 Introduction 

2 Spherical dust collapse 

2.1 Metric, singularity, and apparent horizon 

2.2 Initial conditions for stellar collapse . . . 

2.3 Null geodesic equation 

3 Astrophysical censorship 

3.1 Weak censorship (n = 1 and n = 2) . . . 

3.2 Strong censorship (n = 3) 

4 Conclusion 



1 Introduction 



Singularity theorems have been one of the most important developments in the theory of clas- 
sical general relativity pQ. These theorems ensure that under very general and physically 
reasonable conditions singularities are bound to occur in nature. Such singularities may occur 
in a cosmological model or gravitational collapse of a star. In physical terms a singularity is a 
point of spacetime where curvatures become so large that laws of classical physics do not hold. 
Mathematically it is not a point in the spacetime manifold. To ensure the validity of classical 
laws of physics, a hypothesis that a physically realistic gravitational collapse within classical 
general relativity will never result in a naked singularity, called cosmic censorship [2], was pro- 
posed. By a globally naked singularity we mean the existence of future-directed non-spacelike 
curves from a singularity to a far-away observer, implying the violation of weak censorship. 
If all the future-directed non-spacelike curves fall back into the singularity we only have the 
violation of strong censorship, where the singularity is termed as a locally naked singularity. 

The strong censorship has been shown to be false in the interior of charged black holes in the 
Einstein-Maxwell-scalar system numerically [3] and rigorously [I]. Moreover, there has been a 
growing family of examples where naked singularities form in the spherically symmetric collapse 
of a dust cloud [5j|6] and a perfect fluid with nonvanishing pressure [7] , leaving a subtle problem 
whether these examples are "physically realistic" . See [H|9] for other interesting examples. The 
naked singularities seem to arise not only in various matter models (dust, null dust, perfect 
fluids, scalar fields, Yang-Mills fields) but also under various symmetries (spherical, quasi- 
spherical, cylindrical, numerical studies of axially symmetric configurations). All these models 
exhibit the formation of naked singularities as well as black holes. The numerical studies 
of gravitational collapse have also unveiled intriguing phenomenon of critical behavior [TO] . 
Therefore, the non-existence of a general solution to Einstein equations which led the studies 
on gravitational collapse to specific models has yielded rich dividends. In fact, with so many 
possible counterexamples we cannot take the cosmic censorship as a physical principle but 
only as a rather mathematical question of whether or not these naked- singular spacetimes are 
generic in the space of all possible initial data sets. Together with a general belief that general 
relativity should be replaced by quantum gravity in the Planckian regime, it has to be said that 
the original cosmic censorship is obsolete and should be accordingly modified (e.g., see pTT|[T2] ) . 

If the cosmic censorship is violated in any sense, it would be important whether the violation 
can be tested experimentally or observationally. The possibility of the cosmic censorship vio- 
lation in the particle collider has been recently discussed [13J. Studies of the cosmic censorship 
hypothesis from an astrophysical point of view are becoming an active research subject. In 
particular, many researchers are interested in the relevance and consequence of an overspinning 
Kerr spacetime, which is a naked singularity in classical general relativity [T3] . As for the possi- 
ble formation of naked singularity in astrophysical situations, the result of numerical relativity 
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in its early development suggests that the cosmic censorship is preserved in the collapse of 
the evolved star with large angular momentum under plausible initial conditions. See [H] and 
references therein. Giacomazzo et al. [16] recently studied the fate of a rapidly differentially 
rotating neutron star with the angular momentum greater than the Kerr bound and found that 
generic conditions for such a progenitor do not lead to a naked singularity. 

Naively, in the core collapse of a very massive star or the delayed collapse of a newly born 
neutron star, the behaviour of the collapse is understood as free fall once the dynamical in- 
stability sets in and the collapse proceeds. This suggests that we can approach the collapse 
problem by assuming a pressureless fluid or dust. For this reason, we hereafter will concen- 
trate on the spherically symmetric dust model as a first step. The studies so far have mainly 
concentrated on the issue of local visibility and geodesic analysis near the singularity, although 
there have been a few studies analyzing global structures of naked singularity [HIE]. This in- 
troduces a considerable simplification since the formation of naked singularity depends on the 
matter distribution only in the central region. In fact, one can easily construct both locally and 
globally naked singular models when the central matter distribution admits naked-singularity 
formation, if one is allowed to put the mass distribution in the surrounding region arbitrarily. 

In this paper we analyze the formation of singularity and its causal structure based on the 
astrophysically realistic collapse scenarios in the final stage of stellar evolution, such as the core 
collapse and the delayed collapse. The idea is, with so many examples showing the formation 
of naked singularities, to find a suitable physical condition under which the censorship may be 
preserved and study whether or not the astrophysically realistic collapse satisfies this condition. 
One direction could be to impose additional physical conditions on the models showing possible 
counterexamples. However, the imposition of any new condition other than genericity and 
energy conditions can amount to the restrictions which may be hard to justify. One can say 
that naked singularities are harmless in a sense if they do not disturb an asymptotic observer. 
For the asymptotic observer the violation of strong censorship causes no harm. Therefore, as 
for the direct testability of the violation of the cosmic censorship, it is important whether or not 
the weak cosmic censorship conjecture is preserved or not. With this perspective we analyze 
here the well known Lemaitre-Tolman-Bondi model with the initial conditions where naked 
singularities are known to appear as an outcome of the evolution. 

We begin with reviewing the basic framework of the dust-collapse model. The next section 
deals with the main results of this paper. We work in the geometrized unit (c = G = 1), 
otherwise denoted. 
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2 Spherical dust collapse 



The most general solution describing a spherically symmetric dust ball is given by Lemaitre- 
Tolman-Bondi metric [18]. It is a natural generalization of Oppenheimer-Snyder model which 
describes collapse of a homogeneous dust cloud [T9] . 



2.1 Metric, singularity, and apparent horizon 

The energy-momentum tensor of a pressureless fluid is written as T^ v = pu^u v ', with the energy 
density p and the normalized velocity field (u^u^ = —1). Solving the Einstein equations 
with the spherically symmetric ansatz, the Lemaitre-Tolman-Bondi (LTB) solution is obtained 
in comoving synchronous coordinates (w M = 8%) as 

ds 2 = -dt 2 + - ff^ dr2 + R 2 (t,r)dn 2 , (2.1) 
R 2 = f(r) + (2.2) 

^ = (2 ' 3) 

where X' = d r X, X = d t X, and dVt 2 is the line element of a unit two-sphere. F(r) and f(r) 
are integration constants, and fixed once the initial distributions of mass and velocity of the 
fluid are specified. Their physical interpretation can be seen from equation (12.21) : the left-hand 
side is a kinetic energy; the second term on the right-hand side plays a role of potential energy; 
f(r) is a total energy. Integrating equation (12.21) . we obtain R = R(t,r) in an implicit form 

t-t.(r) = -^G (2.4) 



'F 

where t s (r) is a constant of integration and 



" Arcsin v /i/ y/T 



1/3/2 



< y < 1 



G(y)=<l, y = ■ (2.5) 
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-oo < y < 



(_y)3/2 y 

With the coordinate degrees of freedom available, one can choose initial time at t = and fix 
R(0,r) = r without any loss of generality. Then, the expression for time when a fluid element 
defined by a particular r = const, value plunges into the shell-focusing singularity, which is 
defined by R = 0, is given by 



r 



3/2 / j r 



t s (r) = ^=G ( - J - ) . (2.6) 
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For spherically symmetric systems the condition for the existence of an apparent horizon (i.e., 
an inner boundary of the region containing trapped surfaces), corresponds to the existence of 
two-spheres whose outward normals are null, is given by 

V^Vi? = R = F. (2.7) 

We would like to note here that an apparent horizon depends on the slicing of spacetime. 
However, we are interested in the issue of global visibility, and at the boundary of collapsing 
dust ball the apparent horizon coincides with the event horizon. Thus the apparent horizon, 
which is much easier to locate, should serve our purpose here. 

In LTB models, from equations (12.4ft and ( 12. 7ft . we have the following relation between the 
singularity and the apparent horizon (AH) curve 

tAu(r) = t s (r)-FG(-f). (2.8) 

This equation is crucial in understanding the causal structure of singularity in LTB models. 
Since G is a positive and convex function, the positivity of mass F implies £ah(?") < t s (r) for 
the non-central shells (r > 0). It is only at the centre that the regularity condition on the initial 
data demands the mass function F should be zero and hence the singularity and the apparent 
horizon appear simultaneously, i.e., £ah(0) = t s (0). Thus it is only the central singularity which 
can be a naked shell- focusing singularity. In what follows we shall be interested in the central 
shell-focusing singularity. 



2.2 Initial conditions for stellar collapse 

In this paper, we assume the following simple form of the density profile at the initial time 
t = 0, 




< r < r b 

P(0,r)={ r V \LJ " " , (2.9) 

r > r b 

where p c , L, and r b (< L) are positive constants, and n — 1, 2, or 3. Here, p c is the central 
density; r = r b gives the boundary of the star; L is a parameter controlling the gradient of the 
density as well as n\j\ Since F(r) is twice the Misner-Sharp quasi-local mass, the above ansatz 
of density distribution fixes F(r) as 

i r 

-F(r) = An p(0,r)r 2 dr. (2.10) 



o 



1 Since the density profile of neutron stars has a steep cutoff near the surface as is well known [20j . this simple 
choice of the density profile with a discontinuity at the surface is sufficiently good for the present purpose. 
Moreover, as will be shown later, the trapping of null rays from the singularity by the apparent horizon happens 
enough inside the star. Thus, we argue that the detail of density profile around the surface does not affect our 
conclusion. 
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The LTB spacetime is connected to an outer Schwarzschild spacetime whose mass M is given 
by 

M= l -F{r b ). (2.11) 

For our purpose, it is convenient to choose (M, p c , r&, n) as a set of parameters specifying a 
model. From equations ( 12. 9p . (12.101) . and (12. lip , the remaining parameter L can be given in 
terms of these parameters 

/ 127T V /n 

L - rb {(n + 3)(4n-3M/p c rt)J " (2 ' 12) 

Although the set of parameters (M, p c , r&, n) are almost independent, there are two constraints 
among them for L to be real and satisfy L > r&. From equation (12.121) such constraints are 
found to be 

3M\ 1/3 /n + 3 3M xV3 



In the limit that r b takes the value of left-hand side, L diverges, corresponding to a homogeneous 
density profile. The equality holds when L = r&, corresponding to that the density vanishes at 
r — > r b —. 

In addition to F(r), the LTB solution has another arbitrary function f(r), which is related 
to the initial velocity distribution of fluid elements R(0,r). We assume that all fluid elements 
are at rest -R(0, r) = at the initial time, which is called the momentarily static case. From 
equation fl2.2jl . such a case is realized by taking 

/ = --• (2.14) 

r 

This condition is very different from the usual "marginally bound" models, where collapsing 
shells are at rest at spatial infinity, corresponding to vanishing /. The momentarily static 
condition is astrophysically motivated because the collapse to a black-hole formation will begin 
when a dynamical instability sets in. Such an instability is due to: the electron-capture at a 
static massive stellar core; the photo-dissociation reactions at a core of very massive stars in 
the prompt core-collapse scenario; the cooling and mass fall-back at a newly born neutron star 
in the delayed collapse scenario [2T] . 

2.3 Null geodesic equation 

In order to know the causal structure of spacetimes, it is essential to investigate the behavior of 
geodesies, especially those emanating from the singularity in the present case. Here, we review 
the null geodesic equation in the general LTB spacetime (see, e.g., [22] for full descriptions). 
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We are interested in the null geodesies emanating from the central singularity, which is 
characterized by R = r = 0. In particular, one has to identify the first radial null geodesic 
emanating from the singularity, which is a part of the Cauchy horizon, out of all possible null 
geodesies. The central singularity is a singular point of the null geodesic equation from where an 
infinite family of null rays come out. One needs to parameterize these trajectories appropriately 
to make each curve have a distinct tangent at the singularity. Such a parametrization is known 
to be possible in (R, u)-coordinates rather than (t, r)-coordinates, where u := r a with a constant 
a (> 1) to be fixed for each background so that the null geodesies have finite tangents at the 
singularity (R,u) = (0,0). 

The equation of trajectory R = R{u) for outgoing radial null geodesies is 



dR _ L Vf + A/X \ H(X, 
du \ + / J a 

where a complete list of variables and functions is 



u 



(2.15) 



u 



u 



H(X,u) 
9(r) 



(v 



e 



(1 +p)V2 r 3 ( a - 1 )/ 2 



rF' 

V(r) = — , 



rf 



T) - lp\ X" 2 G{-PX) 
(rj - 3/3/2)G(-p) 

r 3(a-l)/2 ' 

P{r) = y , 



P 



1 

X' 



(2.16) 



P(r) 



pr 



a-l 



A(r) = — . 



It is noted that j3(r) is defined to be zero when / is zero. 

The singularity can be naked if there exists an outgoing null geodesic emanating from it 
with a definite positive tangent X := \im Ut ^ R/u = lim Ut ft_>odR/du. An algebraic equation 
for such a positive tangent is obtained by taking the limit (u, R) — > (0, 0) in equation f)2.15p . 



Xn 



v//o + Ao/X \ H(X ,Q) 



a 



(2.17) 



where the subscript denotes the limiting value at r — > 0. Thus, if there exists a positive root 
X of equation (12.171) for given initial conditions F and /, the null ray emanating from the 
singularity with the tangent Xq is proved to be (a part of) the Cauchy horizon, and therefore, 
the spacetime is locally naked at least. 



3 Astrophysical censorship 

As mentioned above, the nakedness of the central singularity is determined by examining the 
null rays locally around the singularity for given initial conditions F and /. In our models, 
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Table 1: Input parameters (p c , M,r b ,n) and parameters (L, X , R tTlip , R tra , p / Rb) obtained from the 
present analysis, a = 5/3 for n — 1 and a = 7/3 for n = 2. 



Model 


p c [10 15 g/cm 3 ] 


M[M ] 


r b [km] 


n 


L [km] 


X [km 1 -"] 


i?trap [km] 


Rti&p/ Rb 


A 








1 


8.98 


0.213 


0.298 


0.0648 


5.00 


1.50 


7.00 














2 


7.09 


0.0583 


0.00170 


0.000339 










1 


13.3 


0.164 


0.533 


0.0778 


B 


2.00 


2.00 


10.5 














2 


10.6 


0.0341 


0.00379 


0.000502 


C 








1 


16.2 


0.144 


0.556 


0.0656 


1.50 


2.70 


12.8 




















2 


12.9 


0.0263 


0.00320 


0.000346 



introduced in section |2~2| F and / are characterized by the parameters (M, p c ,r b ,n). It can be 
shown that for n = 1 and n = 2 the singularity is naked, while for n = 3 the singularity can be 
either naked or censored, depending on other factors. Thus, we discuss the n = 1, 2 cases and 
n = 3 case separately below. We stress here that the n = 2 profile is most physically plausible, 
since it models the mass distribution of stars in hydrostatic equilibrium. This is because in the 
presence of pressure p = p(p), the pressure gradient can be balanced at the centre with the 
gravitational force, which is proportional to r, only for n = 2. For n = 1, the pressure gradient 
force dominates the gravitational force, while the situation is reversed for n = 3. Therefore, 
we are most interested in the n = 2 case, although we also show the result for the n = 1 and 
n = 3 cases because the LTB spacetime is a solution describing the collapse of a regular dust 
cloud even in these cases. 

3.1 Weak censorship (n = 1 and n = 2) 

Since we have a wide three-dimensional parameter space of (M, p c , r^) only constrained by 
(12 . 131) even if n is fixed, it is practical to take several sets of numerical values of (M, p c ,rb), 
each of which is realized by a typical equation of state of neutron star. We adopt three sets of 
parameters, which we call models A, B, and C (see tabled]). Models A, B and C correspond 
to the marginally stable configurations of neutron stars for an extremely soft, medium and 
extremely hard equations of state, respectively (see, e.g., [20]). 

For a given set of (M, p c ,rb,n), parameter L is determined through equation (I2.12p . The 
power a in equation (I2.16P is fixed so that G(r) has a finite value in the limit r — > 0. The 
results are a = 5/3 for n = 1 and a = 7/3 for n = 2. Solving equation (I2.17P with a obtained, 
a positive root Xq is determined. Then, one is ready to integrate H2 . 1 51) to obtain the trajectory 
of Cauchy horizon R = R(u) with the initial conditions of R\ u =o — and dR/du\ u= o = Xq. 
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Let us focus on model B with n — 2. Qualitative behaviors of the other models are the same 
as this case. The Cauchy horizon numerically obtained is shown in figure [TJ Figures Ufa) and 
DJb) show the Cauchy horizon in the (r, R) and (r, t) planes, respectively, with the apparent 
horizon (12. 8ft . One can see that R increases with r at first along the Cauchy horizon, but 
then encounters the apparent horizon at R = i?trap — 3.79 [m] and eventually plunges into the 
non-central singularity R = at r > 0. In order to see that this trapping point is enough inside 
the collapsing star, we define the proper radius of collapsing star at the time of trapping by 

R b := R(t tTaiP ,r b ), (3.1) 

where t trap is the time of trapping measured in the ^-coordinated In the present case, Rb = 
7.56 [km]. Thus, R tra ,p/Rb — 0.05% and we can say that the trapping happens enough inside 
the star. In other words, the naked singularity cannot be globally naked, preserving the weak 
cosmic censorship hypothesis. A conformal diagram is depicted in figure [2j See table [1] for the 
values of parameters for the other models. 

3.2 Strong censorship (n = 3) 

In contrast to the cases of n = 1 and n — 2, where the central singularity is naked necessarily, 
the nakedness of the singularity for n = 3 depends on the other parameters (M, p c , r b ) (it turns 
out that a = 3 in this case). That is, the algebraic equation (12.171) has no positive root for 
certain parameter regions of (M, p c ,r b ), implying that the collapse results in a black hole. As 
will be shown below, one can see that for any set of (M, p c , r b ) corresponding to an initial 
condition of astrophysical stellar collapse, algebraic equation (12. 1 7f) has no positive root. 

In order to reduce the number of parameters and simplify our arguments, let r b be given by 
the right-hand side of the inequality (I2.13p . which results in L = r b . Furthermore, we assume 
that p c takes the value of model B. Then, the only parameter left to be fixed is the mass M. One 
can check numerically that the algebraic equation ( 12.171) has a positive root for M < O.O688M . 
Integrating equation (12.151) for a mass below this threshold, one can see also that the Cauchy 
horizon is never trapped inside the star (M = O.O688M is a marginal case where the Cauchy 
horizon is trapped just at the surface), implying that the singularity is globally naked. On the 
other hand, for M > O.O688M equation (12.151) has no positive root, implying the formation 
of a black hole. Thus, for the collapse of massive stars whose mass is a few solar mass the 
singularity cannot even be naked, preserving the strong censorship. 

2 Although this proper radius has no coordinate-independent meaning of course, it will do for our present 
purpose. 



s 




r [km] r [km] 

(a) (b) 

Figure 1: (a) The trapping of the Cauchy horizon (solid) by the apparent horizon (dashed) in the 
r-R plane for the model B with n = 2. The Cauchy horizon emanating from the central singularity 
R = r = is trapped by the apparent horizon at (r trap , _R trap ) = (0.673,0.00379), which is enough 
inside the stellar surface located at r = r?, = 10.5, and then plunges into the non-central singularity 
lying on R = at r = 0.790. (b) The trapping of the Cauchy horizon (solid) by the apparent 
horizon (dashed) in the r-t plane again for the model B with n = 2. The shell-focusing singularity 
is drawn by a dotted curve, into which the Cauchy horizon plunges at (r, t) = (0.790,47.052). 



4 Conclusion 

We saw in section 13.11 that the naked singularities appearing in the collapse of massive stars 
cannot be globally naked. Here, we illustrate that globally naked singularities can happen only 
for parameters that are unrealistic from an astrophysical point of view. As in the argument in 
section l3~2l let r& be given by the right-hand side of inequality (12.131) . and we take n = 2 and 
the value of p c of model B. Then, only parameter left to be fixed is the mass M. Integrating 
equation ( 12.151) numerically up to the surface for various values of M, one can easily find that 
for M > O.2466M the Cauchy horizon is trapped by the apparent horizon before reaching the 
surface of star as in the cases in section 13.11 whereas for M < 0.2466M Q the Cauchy horizon 
is not trapped within the stellar interior. Thus, the globally naked singularities can appear 
only for collapse of the small-mass regime that is far below the masses of the observed neutron 
stars (~ 1.5M Q ) and the maximum mass (~ 1.5M ) of the neutron star for the extremely soft 
equation of state. 

In this paper, we have modeled the collapse of massive objects by the Lemaitre-Tolman- 
Bondi solution ( 12. ip with the initial conditions characterized by the initial density profile ( 12. 9p . 
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r>0 i + 




i 



Figure 2: A conformal diagram of the LTB solution mimicking the initial mass distribution of the 
astrophysically realistic collapse in the final stage of stellar evolution. The central (r = 0) singularity 
is locally naked but not globally. Namely, the Cauchy horizon (CH, thick solid) is trapped by the 
apparent horizon (AH, thin dashed) and plunges into the non-central (r > 0) singularity without 
reaching the surface of collapsing star (r = r^). The line EH (thin solid) represents the event 
horizon. 

velocity distribution ( 12. 141) . and the central density, mass, and radius corresponding to the 
several models of marginally stable neutron stars (see tabled]). For n—1 and n = 2 (the latter 
is the most plausible), the singularity in all models is naked necessarily, but we have shown that 
it cannot be globally naked. Namely, the first null ray emanating from the singularity cannot 
reach the future null infinity I + (see figure [2]), implying that the singularity cannot be seen by 
any observers except those who plunge inevitably into the future spacelike singularity. Thus, 
it can be said that the singularity is censored by the astrophysical censor in the weak sense. 
For n = 3, the singularity can be naked or not, depending on the parameters such as mass and 
central density, but we have shown the singularity cannot even be naked for astrophysically 
reasonable parameters. Thus, it can be said that the singularity is censored by the astrophysical 
censor in the strong sense in this case. 

Finally, we would like to stress that we are not in the position to claim that there appear 
locally naked singularities in the astrophysically realistic collapse. The appearance of naked 
singularities in our analysis is entirely because we choose the dust as a matter model and the 
situation can be very different in other choice. Instead, our conclusion is that even if a naked 
singularity could appear at the centre of realistic stellar collapse, it will be hidden behind an 
event horizon and cannot be observed by a distant observer. Since the exposition of the central 
naked singularity to a distant observer depends on the dynamics of not central but surrounding 
region, our conclusion is not very sensitive to the choice of matter models. 
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